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ABSTRACT 

The multi-object fibre-fed spectrograph AAOmega at the Anglo-Australian Tele- 
scope has been used to establish and measure accurate (<1 km s~^) radial velocities 
for a new sample of members in the outer parts of the stellar system uj Centauri. The 
new sample more than doubles the number of known members with precise velocities 
that lie between 25' and 45' from the cluster center. Combining this sample with earlier 
work confirms that the line-of-sight velocity dispersion of uj Cen remains approximately 
constant at ~6.5 km s^^ in the outer parts of the cluster, which contain only a small 
fraction of the total cluster stellar mass. It is argued that the approximately constant ve- 
locity dispersion in the outer regions is most likely a consequence of external influences, 
such as the tidal shock heating that occurs each time cj Cen crosses the Galactic plane. 
There is therefore no requirement to invoke dark matter or non-standard gravitational 
theories. 

Subject headings: globular clusters: general; globular clusters: individual {uj Centauri, 
NGC 5139); stars: kinematics and dynamics 
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1. Introduction 

The stellar system uj Centauri has been known to be unusual, at least as regards its stel- 
lar population, for almost four decades. There is now an extensive body of work which shows 
that, unlike the situation for most globular clusters, the member stars of uj Cen possess a large 
range in heavy element abund ance together with distinctive element-to-iron abundance ratios (e.g., 
Johnson &: Pilachowskil 120101 . and the references therein). A substantial spread in Helium abun- 
dance has also been inferred from the observe d abundance s and structure of th e lower main sequence 
in the cluster color-magnitude diagram (e.g.. lNorrisll2004l : lPiotto et al.ll2005l ). Investigations of the 
metallicities of uj Cen stars in the vicinit y of the main sequeri ce turnoff also suggest that the cluster 
has an age spread of perhaps 2 Gyr (e.g.. |Pancino et al.ll201ll . and the references therein). Together 
these characteristics have led to the suggestion that uj Cen has not evolved in isolation but is instead 
the nuclear remnant of a now disrupted nucleated dwarf galaxy that was accreted by the Milky 



Way (e.g.. lFreemaiilll993l ). iBekki &: Freemanl (|2003l ) have shown that despite the tightly bound and 
retrograde current orbit of uo Cen, such a disruption and accretio n proc ess is dynamicahy plausible. 
Nevertheless, the spectroscopic survey of iDa Costa &: ColemanI (J2008l ) showed that there is little 
evidence for any significant ex tra-tidal populat io n surrounding oj Cen a t the p resent day, consistent 
with the photometric study of lLaw et al.1 (J2003|). IPa Costa fc ColemanI (|2008l . hereafter DC08) give 
an upper limit of 0.7% for the fraction of the cluster mass contained between 1 and 2 cluster tidal 
radii. This result requires the tidal stripping and disruption process of the postulated progenitor 
system to be largely complete at early epoc hs with the stars from the disrupted dwarf galax y now 



widely distributed around the Galaxy (e.g.. IWvlie-de Boer et al.l 120101 : iMajewski et al.ll2012l ) 



While the nucleosynthetic history of uj Cen is complicated and not fully understood, the dy- 
namics of the present-day stellar system, at least for the part of the cluster containing most of the 
stellar mass, ar e relative l y wel l established. There h a ve been a nurnber of models of the system in - 
clud ing those oflMevlanI (IiQStI) . I-Mevlan et all (|l995l l , iMerritt et aP (| 19971 ) , Iciersz fc Heggiel ((2003) 
and Ivan de Marel fc AndersonI (l2010l ). all of which, within their adopte d assumptions, reproduc e 
well the available observational data. The most detailed model is that of Ivan de Ven et al.l (|2006l ). 
This axisymmetric dynamical model, which includes rotation and radially varying anisotropy, sug- 
gests that the mass-to-light ratio of cj Cen does not change with radius - the variation in the model 
MjLy value does not deviate significantly fro m the best-fit constant value of 2.5 (solar units) out 
to the limits of the modelled data at r ~ 20' (jvan de Ven et al.ll2006l ). 



However, it is necessary to keep in mind that models such as that of Ivan de Ven et al.l (J2006l ) 
are constrained by the extent of available observational data. In the case of the velocity dispersion 
profile for lo Cen, the data have been limited, until relatively recently, to a radius of approximately 
20' from the cluster center. While this radius (~4 half-light radii) contains most of the cluster stellar 
mass, it nevertheless is less than half the nominal "tidal radius" of w Cen (57'; see the discussion in 
DC08). The lack of information on the velocity dispersion profile at large radii may mean we are 
currently missing some interesting astrophysics. For example, if cj Cen is the nuclear remnant of 
a disrupted dwarf galaxy then it is possible it has retained some of the dark matter content of the 
original system. One of the best places to constrai n the dark matter content is in the outer parts o f 
the cluster where the stellar densities are low (e.g.. lCarraro &: Liall2000l : iMashchenko &: Sillsll2005l ) . 



Scarpa et al.l (J2003l ) presented the first data for the line-of-sight velocity dispersion of lo Cen 
beyond ~20' from the cluster center. They used accurate radial velocities for 75 members with ~20 
< r' < 30 to show that the cluster velocity dispersion profile may be relatively flat beyond 20'. This 
is in contrast to the monotonically declini ng dispersion p rofile expected for a system in dynamical 
equilibrium in which mass follows light. IScarpa et al.l (120031 ) chose to interpret their results as 
indicating the breakdown of Newtonian dynamics i n a weak accelera t ion re gime. However, this 
interpretation has been questioned by, for example, iBaumgardt et al.l (|2005l ) who argue that the 
external infiuence of the Milky Way on clusters such as a; Cen, that lie relatively close the Galactic 
Center, is sufficiently large that the effective acceleration is larger than the critical MOND constant 
ao; thus Newtonian dynamics should still apply. 



Two additional st udies of the v e locity dispersion in the outskirts of oj Cen have recently 
appeared. In the first, ISollima et alj (|2009l ) conducted a survey for new uj Cen members in the 
outer re gions of the cluster a nd combined their radial velocity results with those from the earlier 
study of iPancino et alJ (|2007l ) to generate a velocity dispersion profile for the cluster that reached 
a radial distance of '^^32^ T he typical uncertainty in the velocity dispersion measures was <1 
km s~^. ISollima et alJ ( 20091 ) claim that the velocity profile decreases monotonically from the 



center outwards thoug h their outermost da ta point lies above the previous point by more than the 
combined (la) errors. ISollima et al.l (|2009l ) note that this occurence mi ght be compatible wi th the 
onset of tidal heating in the outskirts of the cluster. Nevertheless, the ISollima et al.l ( 2009) data 



Sollima et al. 



are not inconsistent with a constant velocity dispersion beyond r ~ 20'. The 

sample contains 98 uj Cen members beyond r ~ 20' but of these stars only 13 lie beyond 30'. 



(120091) 



The se cond recent paper is that of lScarpa &: Falomd (J2010l ) in which the data sets of lSollima et al 



(|2009l ) and IScarpa et al.l (J2003l ) have been combine d to p rovide a further estimate of the velocity 
dispersion profile. The addition of the IScarpa et al.l (|2003l ) velocities increases the number of stars 
with accurate velocities for radial distances between 20' and 30' but does not contribute any new 
members beyond r ?^ 30'. The combined data are consistent with a fiattening of the velocity disper- 
sion beyond r ~ 20'. IScarpa &: Falomd (|2010l ) conclude that this dispersion profile clearly deviates 
from the "Newtonian prediction" , by which they mean a monotonically declining dispersion profile, 
and is "best explained by a breakdown of Newtonian dynamics below a critical acceleration" . 

Clearly there is an urgent need for additional accurate radial velocities for bona-fide members 
of UJ Cen in the outskirts of the cluster, particularly beyond r w 30'. The generation of such a 
sample and a redetermination of the line-of-sight velocity dispersion profile in the outer regions of 
u Cen is the purpose of this paper. The sample selection, the observations, and the measurement 
of the radial velocities are discussed in the next section. Section 3 discusses the membership 
status of the candidates, which is important given the low density of u Cen members in the outer 
regions of the cluster. Section 4 presents the velocity dispersion profiles deriv ed from both the new 
obs ervations and froin comb ining the new data with the velocities given by ISollima et al.l (120091 ) 
and IScarpa &: Falomd (120101 ). The results are presented and discussed in §5. 



2. Observations and Reductions 



2.1. Sample Selection 

The list of candidates to be observed consisted first of the 154 probable uj Cen members 
identified in DC08. These stars lie between 20' and 55' from the cluster center, possess velocities 
and line strengths consistent with cluster membership, and have 15.4 < V < 16.75 (DC08). Then, 
in order to increase the sample of potential members in the outer parts of the cluster, two further 
lists were generated. The first was simply the stars in the original DC08 sample that lie between 
30' and 60' from the cluster center and which were not observed in DC08. There are 545 candidates 



in this category. Note that the outer radius hmit of 60' was chosen to match with the 2 degree 
diameter field-of-view of the fiber-positioning system at the prime focus of the AAT, enabhng fiber 
configurations to be centered on the cluster. It also matches the "tidal radius" of a; Cen. 

The second candidate list was generated from the same photometry set as used in DC08 but 
the selection window was extended '^0.5 mag fainter parallel to the cluster sequence in the color- 
magnitude diagram (see Fig. 1 of DC08). As for the brighter sample, the selected stars lie between 
30' and 60' from the cluster center. There are 1798 additional candidates is this list. 



2.2. Observations 

Five nights in 2008 April were allocated to this program with the 2dF fibre positioner and the 
AAOmega spectrograph on the 3.9m Anglo- Australian Telescope (AAT). The fibre positioner at 
the AAT prime focus can allocate a maximum of 392 fibers within the 2 degree diameter field-of- 
view. Each fibre configuration then consisted of 6-8 guide fibre bundles, ~300 fibres allocated to 
candidate members, up to 50 fibres allocated to blank-sky and ~20 fibres to likely members from 
DC08. These latter stars were incorporated into every configuration to monitor for any systematic 
effects in the radial velocities. The fibres from the positioner are fed to the AAOmega s pectrograph 



a double beam instrument with separate blue and red cameras (jSaunders et al.ll2004l : ISharp et al. 



20061 ). The spectrograph was configured with the A5700A dichroic and the 1500V (blue) and 1700D 
gratings (red). The blue spectra cover the wavelength interval AA4940-5650A at a resolution A/ A A 
of ~4000. The red spectra were centred at A8600A with coverage from A8340A to A8775A including 
the Ca II triplet lines at A8498, 8542 and 8662A. The resolution is ~10,000 and the scale corresponds 
to 8.5 km s~^ per pixel. In this paper we concentrate on the red camera data only. 

Less than ideal weather meant that only 11 a; Cen candidate member fibre configurations were 
observed, principally on 2008 April 25 and 2008 April 26. Two of these configurations were repeats 
to compensate for diminished signal due to cirrus affecting earlier observations. Each configuration 
was observed as a set of 3 x 1500 sec exposures preceded by fibre-flat and arc lamp exposures and 
followed by a second arc lamp exposure. A number of bright radial velocity standards were also 
observed through individual fibres during the run. These provide template spectra for the cross- 
correlation analysis used to determine radial velocities. Each dataset was reduced with the pipeline 
reduction code ^d/dcJ which generates a wavelength-calibrated sky-subtracted spectrum for each 
object fibre. The relative fibre transmissions were set from the data using the SKYLINE(MED) 
option which makes use of the significant flux in each raw spectrum from the numerous bright 
night-sky emission lines in the wavelength region covered by the red camera data. The individual 
spectra from each integration were then median combined to remove cosmic-ray contamination. 

An additional configuration made up of likely new members as determined from the April 



^ www. aao.gov.au/AAO/2df/aaoinega/aaoinega_software. html #drcontrol 



observations was observed on 2008 May 31 via the AAT Service Observing program. The instru- 
mental set-up for the service observations on the red side of AAOmega was the same as for the 
April observations and the data were observed and reduced in an identical fashion. 



2.3. Radial Velocities 

Radial velocities were determined by cross-correlation techniques using the IRAF routine fx- 
cor. The template was a high signal-to-noise spectrum of the y = 7.7 F6V star HD 160043, which 
provides a good match to the w Cen member spectra, particularly as regards the width and depth 
of the Ca II triplet lines. The wavelength interval used for the correlation was AA8470-8740A which 
encompasses the Ca II triplet lines as well as a number of weaker lines, while minimising regions of 
potential significant residual from the sky subtraction. After the cross-correlations were computed, 
the output velocity error and the cross-correlation peak height were plotted against the continuum 
count level in the correlation wavelength region for each of the 12 observed configurations. This 
enabled the identification of occasional situations where the correlation had been affected by in- 
strumental effects such as inadequate cosmic-ray removal. In such cases the problem was corrected, 
usually by interpolating over the effected pixels, and the cross-correlation repeated. In the final 
analysis velocities that had output errors exceeding 5 km s~^ and/or cross-correlation peak heights 
less than 0.5 were discarded - these always coincided with the lowest signal-to-noise spectra. 

After applying appropriate heliocentric corrections, the zero point of the velocity system was 
determined by correlating the template spectrum with other observations of the same standard 
in different fibres and with similar spectra of three other radial velocity standards (HD83516, 
HD101266 and HD162356). The zero point was set by minimising the difference between the 
observed relative velocities for these stars and their catalogue values. The 11 observations of the 
four standards then have a standard deviation of 0.8 km s^^ about their catalogue values indicating 
the velocity zero point is well determined. Only the 245 stars with velocities exceeding 100 km s~^ 
were then retained for the subsequent analysis. 

To determine the velocity errors the rms deviation about the mean velocity was first calculated 
for the ~20 uj Cen members observed in the majority, if not all, of the 12 configurations. Both the 
mean velocity and the rms were calculated using the output velocity errors from fxcor as weights. 
These data show that the rms about the mean velocity is below 1 km s^^ when the continuum 
level in the cross-correlation region exceeds ~900 ADU. The rms values then rise relatively rapidly 
with decreasing continuum levels to 1.7 it 0.3 (la) km s^"*^ at '^600 ADU. A similar analysis of 
the stars observed in the two repeated configurations is consistent with these values and indicates 
further that the rms continues to rise with decreasing continuum level to ~3 it 0.7 (la) km s~^ 
at ~200 ADU, the lowest continuum level of the stars remaining in the data set. For stars with 
a single observation the velocity error was then set by this (rms, continuum level) relation while 
for stars with multiple observations the error was taken as the rms divided by the square-root of 
the number of observations. The overall median velocity error is less than 1 km s~^ excluding any 



systemat ic zero point unce rtainty. This velocity error is lower than the 2-3 km s~^ veloci ty error 
listed bv lLane et alj (|2009l ) who used a similar observing setup, although iLane et alj (|2009l ) do not 
give any information on the continuum levels of their spectra. The uncertainty in the velocity errors 
is sufficiently small that its contribution to the uncertainty in the calculated velocity dispersion is 
negligible. 



3. uj Cen Membership 



In Fig. [T] we plot t he radial velocity, corrected for perspective rotation (see lSollima et al.ll2009l : 
van de Ven et al.l 120061 ) against distance from the cluster center in arcmin. The radial distan ces 



are computed using a tangent plane projection (see ISollima et al. 
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van de Ven et al.ll2006l ) to 



allow for the large angular diameter of the field surveyed. Considering first the 109 stars from 
DC08 with radial distances between 20' and 30', it is evident that the vast majority of these stars 
are apparently probable cluster members despite the relatively low velocity precision (~11 km s~^) 
of the earlier study. Only one star, 8_7_16862, is definitely reclassified as a non-member based on 
the velocity determined here of 278.8 km s~^, while a second star, 8_3_1066 with \j. = 254.2 km 
s~^, lies just outside the ±20 km s~^ from the cluster mean boundaries shown in the figure. This 
star is retained as a possible member for the moment. 

For the 136 stars with v^ > 100 km s~^ and which lie beyond 30' from the cluster center, Fig. [1] 
shows that there is an apparent grouping around the cluster mean velocity out to a radius of at 
least 45' and possibly beyond. However, it must be kept in mind that the surface density profile 
of the cluster is dropping rapidly with increasing radius: the profile given in DC08 indicates that 
the cluster star density drops by a factor of 5 between 20' and 30' and by a further factor of 10 
between 30' and 40'. Conversely, the area that needs to be surveyed goes up as r^, as does the 
number of contaminating non-members assuming they have a uniform surface density. For these 
reasons coincidence with the uj Cen mean velocity does not guarantee cluster membership in the 
outer regions of the cluster, and additional information must be used to help exclude non-members. 
Fortunately, we can make use of the known properties of the stellar population of u Cen to carry 
out this task. 



The stellar system oj Cen is well known for its internal spread in [Fe/H] abundance. iJohnson &: Pilachowski 



(J2010l ) have provided [Fe/H] values for a large sample of a; Cen red giants and their results show, 
in agreement with earlier work (see lJohnson &: Pilachowskill2010l . for references), that: (a) there is 
a lower bound at [Fe/H] ~ -1.9 to the abundances of w Cen red giants that remains constant with 
increasing radius; (b) there is a decrease in the number of stars with [Fe/H] > -1.3 with increasing 
radius relative to the number of more metal-poor stars; and (c) t he dispersion in abundance fo r 
stars with [Fe/H] < -1.3 is approximately constant with radius. The I Johnson &: Pilachowskil (J20ld ) 
sa mple reaches out to radial distances of only ~24' but the limited sample of more distant stars 



m 



Norris et al.l () 19971 ) suggest that these results apply also at larger radii. Further, despite the 



complexity of the oj Cen color-magnitude diagram (CMD) in the vicinity of the main sequence 




Distance from Cluster Centre (arcmin) 

Fig. 1. — The radial velocity, corrected for perspective rotation, is plotted against distance from 
the center of uj Cen for all stars observed whose velocities exceed 100 km s~^. Stars between 20' 



and 30' come from the sample of probable members of lDa Costa &: ColemanI (J2008l ) while the stars 
at larger radii includes stars from that sample as well as stars observed here for the first time. The 
dashed line is at the cluster mean velocity of 233 km s~^ and the dotted lines are at ±20 km s~^ 
from the mean. 



turnofF (e.g. iBellini et al.ll20ld ). it is likely that the age-r ange among the cluster stars is compara- 
tively small, less than 2-4 Gyr (e.g.. iPancino et al.ll201ll . and the references therein). Thus there 
should be a reasonable level of consistency between overall abundance inferred from location on 
the red giant branch in the CMD and that inferred spectroscopically. We can therefore use the 
photometry and the line strengths for the stars in the velocity window 213-253 km s~^ to select 
probable members at all radial distances. 

The approach is as follows. In Fig. [2] we show in the upper panel the CMD for the 108 stars 
in the velocity range 213-253 km s~^ (plus 8_3_1066) which have radial distances between 20' and 
30'. The lower panel shows the CMD for stars in the same velocity range but with radial distances 
between 30' and 60'. These two groupings are designated the inner and outer samples, respectively. 
Note that for completeness we have included in the outer sample stars 5_3_226 and 9_4_1918 as with 
velocities of 212.5 ± 1.2 and 211.4 it 1.5 km s~^, respectively, they lie very close to the lower limit 
of the velocity selection range (see Fig. [1]). Show n also in both pane ls are theoretical isochrones 
from the Dartmouth Stellar Evolution Database (jDotter et al.l |2008| ) for an age of 13 Gyr and 
metallicities [Fe/H] of -2.0, -1.5, -1.0 and -0.5 dex, respectively. For the three lower metallicities 
the isochrones are for [a/Fe] = +0.4 while the most metal-rich isochrone has [a/Fe] = + 0.2 dex. 



Pancino et al. 



2002 



This is consistent with the de pendence of [a/Fe] on [Fe/H] in the cluster (e.g 

Johnson fc Pilachowskil |2010| ). The isochrones have been fitted assuming (m-M)y = 13.94 and 

E{V - I) = 0.16 mag (E(B - V) = 0.1 2 mag) as tabulated in the 2010 version of the Milky Way 



Globular Cluster database (JHarrislll996l . hereafter HIO). 

Similarly we show in Fig. [3] plots of the combined equivalent widths of the A8542A and 
A8662A lines of the Ca II triplet against V — Vhb and V — I for the inner and outer samples. 
The equivalent widths have been determined using gaussian fits to the line profiles with feature 
and co ntinuum bandpasses similar to those adopted in the original work of lArmandroff &: Da Costa 
(Il99ll ). The Vhb value for oj Gen was taken from HIO. The values oi V — Vhb he outside the 
range of existing abundance calibrations for the Ca II triplet which are tailored to more luminous 
red giants. Nevertheless we can use the inner sample, which is dominated by cluster members, to 
define reasonable upper and lower envelopes for the relation between line strength and V — Vhb 
followed by uj Gen stars. The adopted linear relations are shown as the dashed lines in the upper 
left panel of Fig. [3l The relations are then duplicated for the outer sample as shown in the upper 
right panel of the figure. In both cases the lower envelope is very well defined consistent with the 
result for more lum inous samples that the lower abundance cutoff i n the abundance distribution 
is quite sharp (e.g., iNorris et al.lll996l : I Johnson &: Pilachowskil |2010| ). The upper envelope is less 
well define d, as expected, gi ven the significant range of metallicities present even at large radial 
distances (JNorris et al.l 119971 ). The lower panels of Fig. [3] show the relation between the Ca II 
triplet strength and {V — I) colour for the two samples. Members of oj Gen should show a broad 
correlation between these quantities with, in general, stronger lines going with redder colours. The 
panels show this is the case. 



We then combine the information in Figs. [2] and [3] as follows. For each star in the panels of 
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Fig. 2. — Color-Magnitude Diagrams for stars with radial velocities between 213 and 253 km s~^ 
and distances from the cluster center of 20'-30' (upper panel) and 30'-60' (lower panel). Shown 
also are Dartmouth isochrones for an age of 13 Gyr and (left-to-right) ([Fe/H], [a/Fe]) values of 
(-2.0, +0.4), (-1.5, +0.4), (-1.0, +0.4), and (-0.5, +0.2), respectively. The isochrones have been 
fitted assuming (m-M)y = 13.94 and E(y — I) = 0.16 mag. The stars adopted as cluster members 
are plotted as filled symbols while the likely non-members are plotted as open symbols. 
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Fig. 3. — The sum of the equivalent widths in A of the Ca II triplet lines at A8542A and A8662A lines 
is plotted against V — Vhb (upper panels) and V — I (lower panels) for the inner sample (left panels) 
and the outer sample (right panels). The linear segments in the upper left panel encompass the 
range of line strengths shown by probable cluster members in the inner sample. The segments are 
replicated in the upper right panel. The stars adopted as cluster members are plotted as filled 
symbols while the likely non-members are plotted as open symbols. 
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Fig. [2] we measure the offset between the V — I color of the star and the color of the 13 Gyr, [Fe/H] 
= -2.0, [a/Fe] = +0.4 isochrone at the star's V magnitude, normalized by the color difference 
between this isochrone and that for 13 Gyr, [Fe/H] = -0.5, [a/Fe] = +0.2 at the V mag of the star. 
For the same star we then measure the difference in equivalent width between the value for the star 
and the lower envelope line shown in the upper panels of Fig. [3] at the V — Vhb value for the star. 
This difference is then normalized by the equivalent width difference between the upper and lower 
lines at the V — Vhb value of the star. For members of the cluster these two quantities, denoted 
by ^^(SW) and 5n{y — I) respectively, should be well correlated as both are metallicity indicators. 
Non-members, however, will in general lie away from the cluster member sequence. The results of 
this process are shown in Fig. U] where the expected correlations are evident. The dashed lines in 
the upper panel of the figure show the membership selection window adopted for the inner sample. 
The selection window is then reproduced in the lower panel for the outer sample. 

We have then combined the information from Figs. [2l [3] and [5] to provide our best estimate of 
the uj Cen membership status for the stars in the inner and outer samples. The adopted cluster 
members are plotted as filled symbols in all three figures while the non-members are plotted as open 
symbols. We note that the membership status of stars falling near the boundaries of the selection 
window in Fig. [J] were individually considered and classified taking into consideration uncertainties 
in the photometry and line strength measurements. In particular, despite our efforts to minimise 
them it is still likely that there are systematic uncertainties in the {V — I) photometry at the ±0.03 
mag level. For inner sample, 93 of the original 108 stars are classified as members, including star 
8_3_1066, while for the outer sample, 67 of the original 91 are classified as members, including 
5_3_226 but not 9_4_1918. We then list in Table [U the identification, J2000 position, heliocentric 
velocity and error, distance from the cluster center in arcmin, V and V — I photometry, and the 
sum of the equivalent widths of the A8542A and A8662A Ca II triplet lines in A, together with 
its associated error, for the 160 adopted cluster members. Table [2] gives the same information for 
39 stars from the inner and outer samples that are classified as probable non-members of uj Cen. 
Figure [5] shows the observed velocity versus radial distance diagram of Fig. [1] but now with the 
probable members and non-members identified. 



We have then used the Besancon model of the Galaxy (JRobin et al.l 120031 ) to estimate the 
success of this membership selection process. Five independent realisations of the model Galaxy 
were generated for the line-of-sight towards uj Cen using an area on the sky equivalent to that of 
the 2dF field-of-view. The V and V — I magnitude and color ranges for the model output were 
chosen to match approximately those of the outer sample of stars (see the lower panel of Fig. [2]). 
Normalisation of the model to the observational data was then set by the ratio of the number of 
stars in the outer sample with velocities in the range 140-190 km s"^ to the number of model stars 
in the same velocity interval. The predicted number of field stars in the velocity interval 213-253 
km s~^ could then be calculated from the model numbers. The predicted number of field stars with 
velocities exceeding 260 km s~^ was also calculated as a check. We find that the model normalised 
in this way over predicts the number of high velocity stars: 8+1 stars are predicted versus the 4 
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Fig. 4. — The normalized relative line strength (5„SW in A is plotted against normalized relative 
giant branch color dniV — I) for the inner sample (upper panel) and for the outer sample (lower 
panel). The dashed lines in the upper panel outline the membership selection criteria adopted 
for the inner sample. The lines are reproduced in the lower panel. The stars adopted as cluster 
members are plotted as filled symbols while the likely non-members are plotted as open symbols. 
Stars near the boundaries were considered individually and classified taking into consideration 
uncertainties in the photometry and line strength measurements. 
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Table 1. uj Cen Probable Member Data 



ID 


RA (2000) 


Dec (2000) 


V. 


^(V,) 


r' 


V 


V -I 


sw 


e 








(km s^i) 


(km s-^) 




(mag) 


(mag) 


(A) 


(A) 


8_7_15831 


13 25 06.54 


-47 17 41.6 


223.4 


0.4 


20.0 


15.80 


1.01 


2.99 


0.09 


8_4_3206 


13 27 47.33 


-47 45 44.2 


241.6 


1.0 


20.1 


16.42 


0.97 


2.25 


0.10 


8_8_3134 


13 25 47.75 


-47 10 59.8 


232.6 


0.7 


20.2 


16.69 


1.02 


2.61 


0.11 


8_6_16385 


13 24 56.26 


-47 36 34.7 


235.7 


0.6 


20.2 


16.38 


1.02 


1.95 


0.07 


8_2_1336 


13 28 43.30 


-47 24 37.8 


234.8 


0.4 


20.2 


16.19 


1.04 


2.58 


0.09 



Note. — This table is available in its entirety in a machine-readable form in the online journal. A 
portion is shown here for guidance regarding its form and content. For stars in common these data 
supercede those in Table 2 of Da Costa & Coleman (2008). 



Table 2. uj Cen Probable Non-Member Data 



ID 


RA (2000) 


Dec (2000) 


V, 


^(V,) 


r' 


V 


V -I 


sw 


e 








(km s-i) 


(km s'l) 




(mag) 


(mag) 


(A) 


(A) 


8_8_2219 


13 26 04.35 


-47 09 39.7 


224.8 


0.7 


20.2 


16.40 


1.09 


3.95 


0.13 


8_8_4146 


13 25 24.37 


-47 13 35.6 


226.8 


0.8 


20.4 


16.15 


1.12 


3.05 


0.10 


8_1_2780 


13 27 51.11 


-47 11 28.5 


235.7 


0.3 


20.4 


15.83 


1.02 


3.85 


0.18 


8_8_4052 


13 25 26.75 


-47 12 55.7 


240.5 


0.4 


20.6 


16.13 


1.08 


2.15 


0.07 


7_3_412 


13 24 42.89 


-47 30 51.8 


244.6 


0.4 


20.9 


15.99 


0.90 


2.06 


0.07 



Note. — This table is available in its entirety in a machine-readable form in the online journal. 
A portion is shown here for guidance regarding its form and content. For stars in common these 
data supercede those in Table 2 of Da Costa & Coleman (2008). 
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Fig. 5. — As for Fig. [T] the radial velocity, corrected for perspective rotation, is plotted against 
distance from the center of lo Cen for all stars observed whose velocities exceed 100 km s^^. 
Adopted probable members are now plotted as filled symbols while probable non-members are 
plotted as open symbols. 



15 



actually observed in the outer sample. The difference however, is not very significant given the 
small numbers. 

For the cluster velocity range, the normalised model predictions are 3=bl field stars for distances 
from the cluster center between 30' and 40', 2±1 field stars for between 40' and 50', and 3±1 field 
stars for between 50' and 1°. The number of stars classified as probable non-members in the 
corresponding radial distance intervals (see Fig. [5]) are 13, 6 and 4 stars, respectively. Assuming 
the validity of the model and the normalisation, the comparison then appears to indicate that the 
membership classification process adopted has in all likelihood been too conservative: perhaps as 
many as dozen of the "probable non-members" are in fact likely to be cluster members, with most 
of the mis-classified stars falling in the 30'-40' radial range. The effect of this potentially overly 
conservative membership selection on the calculated velocity dispersions will be discussed in the 
following section. 

In Fig. [1] there are eight stars with velocities between 213 and 253 km s^^, plus stars 5_3_226 
and 9_4_1918 which lie just outside the velocity interval, that have distances from the cluster center 
exceeding 46'. Only five of these ten stars survive the cluster membership analysis and are listed in 
Table[T]and shown as filled symbols in Fig.[5j This number is too small for a statistically meaningful 
measure of the velocity dispersion at this extreme outer region and so the subsequent analysis will 
be based on the 155 probable oj Cen members that lie between 20' and 46' from the cluster center. 

Before discussing the velocity dispersion profile defined by these data we show in Fig. [6] the 
surface density profile for cj Cen. The data are taken from DC08 (see references therein) except 
that we have used the surface density values determined from the cluster member sample derived 
here in place of the equivalent data in DC08. The surface density points from the inner and outer 
samples have been separately scaled vertically to match the existing data. The separate scaling 
is necessary since the outer sample covers a larger V magnitude range than the inner sample. It 
is worth noting that the outer parts of this profile do not obviously exhibit the increased profile 
steepening at large radii characteristic of a tidally limited profile. Instead it appears that the outer 
profile is best described by a constant slope, i.e., a power-law profile. A (unweighted) least-squares 
fit to the data points lying beyond 20' (1.3 in log r) yields a value of -5.4 it 0.2 for the power-law 
slope. The fit is shown as the straight line in the figure. It is also worth noting that the surface 
density implied by the 5 possible members beyond 45', whose radial distances range between 47.4' 
and 54.0', is consistent with the observed profile and the power-law fit. 



4. The Velocity Dispersion Profile 



The outer sample of ISollima et al.l (|2009l ) contains 98 uj Cen members that have distances 



from the cluster c e ntre e xceeding ^20', althou g h the re are only 13 members beyond 30'. Similarly. 



Scarpa fc: Falomd (120101 . see also IScarpa et al.l (J2003l )) give radial velocities for 75 oj Cen members 



with radial distances beyond 20' from the cluster center, but there is only one member (just) 
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Fig. 6. — The (circula rly averaged) surface density profile for u Cen taken from 



Da Costa &: ColemanI (J2008l ) except that the blue filled circles, which represent the cluster member 
sample determined here, replace the equi valent data in the ear l ier w ork. The vertical dot-dash 
line indicates the tidal radius adopted by iDa Costa &: ColemanI (J2008l ) . The dashed line and the 
vertical arrow indicate the up per limit on the densit y of cl uster members for the region between 
1 and 2 tidal radii derived by iDa Costa &: ColemanI (J2008l ) , while the dotted lines represent the 
statistical uncertainty in the upper limit. The red solid line is a power law fit to the outer points; 
it has a slope of -5.4 dex/dex. The fit to the outer points is enlarged in the insert box. 
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beyond 30^ in that sam ple. There are 26 stars in common between the ISoUima et al.1 (|2009l ) and 
Scarpa &: Falomd (J2010|) samples. The current sample, with 62 probable members between 30' and 
46' therefore represents a considerable increase in the number of velocities with which to study the 
velocity dispersion profile in the extreme outer parts of this stellar system. 



4.1. A AT sample 

In Fig. [7] we show the velocity dispersion points calculated from the present sample of u Cen 
probable memb ers. The velocity disp ersions have been calculated using a maximum likelihood 
estimator (e.g.. iPrvor &: Mevlanlll993l ). In making the calculation the stars have been grouped 
into bins containing at least 10 members - the radial range of each bin and the number of stars 
included are indicated in the lower part of the figure. The dispersions have been calculated relative 
to the mean velocity for each group - these individual mean velocities differ by at most 2.3 km s~^ 
from the mean (233.4 km s~^) for the 93 members with distances between 20' and 30' from the 
cluster center. Calculating the dispersions for the stars beyond 30' relative to this fixed mean causes 
only a slight increase (0.1 to 0.4 km s~^) in the dispersions but the differences are well within the 
uncertainties. 

We also note that in the previous section it was suggested that the membership selection 
employed had been too conservative in that perhaps as many as a dozen of the "probable non- 
members" may be actual cluster members. To investigate the effect of this possibility on the 
calculated dispersions, we conducted five trials in which 12 stars were randomly selected from 
the set of 19 probable non-members that lie in the radial range 30' to 46'. The selected stars 
were then combined with the probable members in the appropriate radial bins and the dispersions 
recalculated. In all cases the change in the dispersion was less than the errors calculated for the 
probable members only samples. The mean change was an increase in dispersion of 0.3 km s~^ with 
the largest excursions seen being an increase in the dispersion in the 33' -36' bin of 1.2 km s~^ and 
a decrease of 0.6 km s~^ in the dispersion for the 40'-46' bin. Both these changes are within the 
error for the equivalent probable members only sample. We conclude therefore that the dispersion 
measurements are stable against modest changes in the membership status of individual stars. 



Shown also in Fig. [7] are the velocity dispersion meas urements from lSollima et al.l ( 20091 ). taken 
directly from their Table 1, fromlvaii de Ven et al.l (120061 ) as the line-of-sight dispersions from their 
Figure 8, and from lScarpa &: Falomd (J20ld ). In the latter case the d ispersion points were calc ulated 
directly from the heliocentric radial velocities listed in Table 2 of IScarpa &: Falomd (J20ld ) using 
identical techniques, including correction for perspective rotation, as those employed for the uj Cen 
stars observed here. 

It is apparent from Fig. [7] that within the radial range where the different samples overlap, 
the velocity dispersion measures are consistent with one another. It is also evident that there is no 
evidence for any significant decline in the line-of-sight velocity dispersion of w Cen members beyond 
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~25'. In particular, the new data, which extend well beyond the previous data, are consistent with 
a constant line-of-sight velocity dispersion of '^6.5 km s~^ in the outer parts of the cluster. 



4.2. Combined Sample 



The 26 stars in common between th e ISollima et alJ (l20C)9l'l an d the IScaroa &; Falomd (J2010|) 



Sollima et al.1 (2009|)) of 0.2 



3)) 



samples have a mean velocity difference ([Scarpa &: Falomd (J20ld ) 
km s~^, with a standard deviation of 0.8 km s~^, confirming the velocity precision of both data sets. 
Unfortunately, due to the difference in the appare nt magnitudes of the samples, there are no stars 
in common between the present work and that oflSollima et al.l (120091'). There are, however, three 



stars in the present sample that are also in [Scarpa fc Falomd (J2010l ). [Scarpa &: Falomd ([2010l ) stars 
0006, 78004 and 85007 corres pond to stars 8,8 4776, 8,5-6453 and 7_4_160 here. For the first two 
stars the velocity differences ([Scarpa &: Falomd ([2010l ) - present work) are gratify ingly small: -0.6 
and -0.2 km s~^, respectively, but for the third star the [Scarpa &: Falomd ([2010l ) velocity exceeds 
that found here by 11.7 km s~^. Star 7_4_160 is one of the stars observed in all 12 configurations 
here and there is no indication of any velocity variab ility: the rms about the w eighted mean velocity 
is 1.2 km s~^. Consequently, we have not used the [Scarpa &: Falomd ( 20ld ) velocity for this star. 
Otherwise, the velocities for the stars in common have been averaged. The combined data set then 
has a total of 299 stars, 224 with radial distances between 20' and 30' and 75 lying between 30' and 
46'. 

Figure [8| then shows the velocity dispersion profile for the combined sample. Shown also 



i n the figure are velocity dispersion measurements calculated from the sample of ISollima et al 



(|2009l ) for stars betwee n 12' and 20'. Th e data are tabulated in Table [31 The line-of-sight velocity 



dispersion points from 



van 



de Ven et al.l (l2006[ ) shown in Fig. [7| are also reproduced in the figure. 
The combined sample clearly verifies what was already evident from Fig. [7| - that although the 
velocity dispersion decr eases outwards with increasing radi us for the inner pa r ts of the cluster (see, 
for example, figure 8 of [Sollima et al.[ ([2009[ ) or figure 4 of [Scarpa Sz Falomd ([20ld )). beyond ~25' 
the velocity dispersion profile shows no signs of decreasing with increasing radius. This is despite, 
as noted above, the surface density dropping by a factor of --^10 between radial distances of 30' 
and 40'. For the 140 stars in the combined sample with r > 25', the mean radius is 31.3' and the 
velocity dispersion is 6.6 ± 0.4 km s~^. 



Discussion 



The first question to be addressed is whether there is a dynamical model which can reproduce 
the surface brightness/surface density profile of the cluster, and the observed line-of-sight velocity 
dispersion profile, without any requirement for dark matter (i.e., a model in which m ass follows light ) 
or similarly, without any requirement for non- Newtonian gravity. For example, the [WilsorJ ([l975[ )- 
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Fig. 7. — The line-of-sight velocity dispersion of u Cen member stars is plotted against distance 
from the cluster center in arc minutes. The filled symbols are for the current probable members 
data set. The radial range corresponding to each dispersion point is indicated in the lower part 
of the plot as is the number of st ars in each radius bin. Shown also on the plot are the l ine-of - 
sight velocity dispersion poin ts of Ivan de Ven et al.l (120061) (open diam onds) , IScarpa et al.l (J2003l ) 
(from lScarpa &: Falomdl2010l . open triangles), and lSollima et al.l (J2009l ) (open circles). There is no 
evidence for any significant decline in the velocity dispersion beyond ~25'. 
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Table 3. Velocity dispersion profile for the outer parts of oj Cen from the combined sample 



Radius Range 


Mean Radius 


N 


0" 


error 


(arcmin) 


(arcmin) 




(km s-i) 


(km s"^) 


12-14 


12.88 


42 


9.62 


1.06 


14-16 


15.12 


44 


8.85 


0.95 


16-18 


16.87 


61 


7.48 


0.69 


18-20 


18.98 


39 


8.49 


0.97 


20-22 


20.98 


89 


7.22 


0.55 


22-24 


23.08 


55 


7.92 


0.77 


24-27 


25.28 


51 


6.10 


0.61 


27-30 


28.55 


29 


7.74 


1.03 


30-33 


31.15 


32 


6.79 


0.87 


33-36 


34.14 


15 


6.03 


1.12 


36-40 


38.03 


16 


6.52 


1.18 


40-46 


43.35 


12 


6.22 


1.32 



Note. — The first four entries are based on the data set of 



Sollima et al.l (|2009l ) while the remaini ng entries are drawn 



fron i the combined data set o f this work, ISollima et al.l (|2009l ) 
and IScarpa fc Falomol (|20ld ). 
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Fig. 8. — The line-of-sight velocity dispersion of a; Cen member stars is plotted against distance 
from the cluster center in arc minutes. The filled circles are for the combined data set of this 



of Sollima et al 



w ork. ISoUima et al.l (120091 ) , and IScarpa &: Falomd (|2010l ) . The filled stars are from the data set 



(I2OO9II for stars between 12' and 20'. Line-of-sight velocity d ispersion points froin 



van de Ven et al.l (|2006l ) are shown as open diamonds. For the combined and ISoUima et al.l (|2009l ) 
points the radial range corresponding to each dispersion point is indicated in the lower part of 
the plot, as is the number of stars in each radius bin. As in Fig. [71 there is no evidence for any 
significant decline in the velocity dispersion beyond ~25'. 



22 



type model for uj Cen presented in lSollima et al.1 (|2009l ) fits the projected surf ace density pro f ile, th e 
ellipticity profile and the rotation curve adeq uately, as shown in Figure 9 of ISollima et alj (|2009l ). 
However, while the model also reproduces the lSollima et al.1 (J2009l ) velocity dispersion data, it does 
not fit the more extensive velocity dispersion profile data presented here. The velocity dispersion 
profile of the model is consistently below the observed points in Fig. [8] beyond ~25', and declines 
monotonically to, for example, 3.6 km s~^ at ~45', significantly below the observations. The model 
is therefore not an adequate description of the dynamics in the outer parts of the cluster. 



Th is is likely to be the case for all s imilar models, e.g.. lMevlan et al.l(|l995l ): lMcLaughlin &: Mevlan 

( 20031 ): iMcLaughlin Sz van der Marell ( 2005l rl. primarily because they are based on a funda mental 



assumption that t he velocity distribution function f(v) is of a 'lowered Maxwellian' form (JKing 



1966 



Wilsonlll975l ). This ensures that the density reaches zero at a fi nite radius, usually identified 



with the boundary set by the tidal force of the Milky Way (lKing]ll966l ) . In such models the velocity 
di spersion profile also declines monotonically, reaching zero at the same finite radius (e.g.. Fig. 1 
of iMeylan et al.lll995l ). For oj Cen, however, the lack of an obvious tidal radius cutoff signature 
in the surface density profile (see Fig. [6]) suggests that models of this type are not appropriate for 
the outer parts of the stellar system. In this sense the disagreement with the velocity dispersion 
observations reveals the inadequacy of the models, not necessarily anything more fundamental. 



A more heuristic approach to the modelling of cj Cen is that taken by Ivan de Mar el &: Anderson 



(|2010l ). whose primary aim was to place constraints on the possible presence of an intermediate- 
mass black hole at the center of the cluster. Their approach was to parameterise the surface 
density profile and then solve the spherical anisotropic Jeans equation to predict the velocity dis- 
persion profiles. The surface bri ghtness profile was fit with a so-called "generalised nuker" profile 



van de Marel &: Andersonll2010l ). which allows for a central power-la w cusp and which contains 



two c haracteristic logarithmic slopes with associated "break" radii (see Ivan de Marel &: Anderson 
2010l . for details). This adopted functional form does not have any particular physical significance, 
but i n this conte xt it is i mpor tant to note that a density profile of this type, in contrast to the 
King (j 19661 ) and IWilsonI (jl975l )-type models described above, does not have a finite cutoff radius 
at which the density, and the velocity dispersion, go to zero. 



W e show in the left panel of Fig. [9]the best-fit "generalised nuker" model of lvan de Marel &: Anderson 
(2010|) for the case where the central logarithmic slope 7 is fixed at zerqfl, compared to the surface 
density data presented in Fig. [6l The model data have been scaled vertically to correspond to the 



^The latter two references show a observed velocity dispersion point of 3.5 ± 1.5 km s~^ at a radius o f 36.6', 
which is in accord with the predictions of the Wilson-type models. This datum conies from ISeitzen i 19831 ) and is 
based on velocities for fou r stars whose r adial distances extend between ~3l' and ~44'. Three of the four stars are 
included in the sample of ISoUima et all ( 20091 ) . We assert that the considerably more extensive data set presented 
here supercedes this early result. 

^This model differs from the overall best-fit model, which has a shallow cusp with 7 — 0.05, only in the very inner 
regions of the cluster (r < 20"). 
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adopted surface density scale of the observations. The fit is excehent at all radii: in particular, the 
relatively constant logarithmic slope at large radii is in agreement with the observations, noting 
that t he new surface density da t a pre sented here, and those given in DC08, were not included 
in the Ivan de Marel &: Anderson! ((20101) fitting process. The right panel shows the corresponding 
line-of-sight velocity dispersion profile compared with the observations from Fig. [HI The velocity 
distribution is mildly anisotropic with the transitio n from radial anisotropy at smal l radii to tan- 
gential anisotropy at large radii occu r ring a t r « 12' ()van de Marel &: Andersonll2010l ) in agreement 
with the results of Ivan de Ven et al.l (|2006l ). The model velocity dispersion curve is consistent with 
the observations over the entire radial range depicted , including the points beyo nd ~25'. For com- 



pleteness we note that the upper panel of Figure 7 of Ivan de Marel fc AndersonI (J20ld ) shows that 
this model also reproduces satisfactorily the observed velocity dispersion profile in the inner parts 
of the cluster. The model explicitly assumes 'mass follows light' and yield s a l^-band mass-to-light 
ratio of 2.6 in solar units, in agreement with the M/Ly value found in the Ivan de Ven et al.l (|2006l ) 
study. These dynamical values are entirely consistent with the mass-t o-light ratio e xpected, given 



plausible assumptions about the stellar population of the cluster (e.g.. lMevlanlll987l ). We can con- 
clude therefore, based on this parameterized model, that an extended dark matter distribution is 
not required to reproduce the observed velocity dispersion profile in the outer parts of the cluster. 
Similarly, given that the modelling process is based on standard Newtonian dynamics through the 
use of the Jeans equation, no non-standard dynamics are required. 



The Ivan de Marel &: AndersonI (|2010l ) model may well be a reasonable description of the dy- 
namics in the outer parts of u Cen but because of its heuristic nature it does not provide any direct 
insight into the physical processes responsible for the applicability of the model. To investigate 
this we note again that the outer surface density profile of u Cen is well represented by a power- 
law and does not exhibit the characteristic "King-profile " tidal cutoff see n in many cluster s (e.g. 



McLaug 



ilin &: van der 



see also 



Marelll200,4I.Tordi fc Grebellboiol ). As discussed bv lPeharrubia et al.l (|2009l . 



Kiipper et al.l ( 2010l )) this is an indication that the phase space in the outer parts of the 
cluster is likely populated up to -E '^ 0, requiring a source of additional energy. The likely source 
of the required heating is the tidal shocks that occur each time uj Cen crosses the Galactic plane, 
as well as the tidal heating the cluster experiences as it moves in the spatially varying potential of 
the Galaxy. 



Dinescu et al.l (119991 ) used the current position and motion of uj Cen to characterise the or- 
bit of the cluster around the Galactic Center. They found that the system has peri- and apo- 
Galactocentric dista nces of approximately 1.2 and 6. 2 kpc, and an orbital p eriod of ~120 Myr 



([Dinescu et al.l [l999|) . Using these orbital parameters Ivan de Ven et al.l ([20061 ) calculate that the 
velocity component perpendicular to the Galactic plane v± is of order 40 km s~^. Consequently, 
for a dis k scale height of 7"250 pc, it takes co Cen about 12 Myr to cross through the disk o f the 



Galaxy ([van de Ven et al.ll2006l ). In contrast, in the to Cen model of Ivan de Ven et al.l ([20061 ) the 
orbital timescale for member stars in the outer parts of the cluster is approxii nately 100 Myr at r 



25' - 30', and longer at larger radii. Thus the impulse approximation (e.g., iGnedin et al. 



199S 



24 





og r(arcrnln) 

Fig. 9. — Left panel: The surface density profile for uj Cen, froi n Fig. [6l The red solid line is the 
surface density profile of the best-fit "generalised nuker" model of lvan de Marel &: Anderson! (J20ld ) 
for the case with the central logarithmic slope 7 fixed at zero. The model profile has been scaled 
vertically to fit the surface density observations. Right panel: The line-of-sight velocity dispersion 
of oj Cen for the outer parts of the cluster. Symbols are as for Fig. [8l The red solid line is the line- 
of-sight velocity dispersion profile of the same "generalised nuker" model whose surface brightness 
profile is shown in the left panel. 
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Binnev &: Tremaindll987l . p. 446) is valid for calculating the "shock heating" the outer parts of the 
cluster experience each disk crossing. 



W e use equation 7-71 of lBinnev &: Tremaind (|l987l ). with parameter values from Ivan de Ven et al 
( 20061 ) ■ to show that the impulsive change in the velocities of stars, \Av\, is '^ 0.36 r' km s~^. This 



change is then comparable to, or exceeds, the line-of-sight velocity dispersion for radial distances 
beyond r ~ 20'. The relative importance of disk shock heating is then measured by a compari- 
son of the shock heating t imescale tghock with the dynamic al timescale tdyn. ■ Use eq uation 7-72 of 
Binnev fc Tremaind (jl987l ) with the parameter values from Ivan de Ven et al.l (|2006l ) gives: 



tshock = iJba^r^ (1) 

for tshock in Myr, a in km s~^ and r in arcmiiu- Adopting a = 6.5 km s~^ for r > 25' and 
then c omparing the radial variation of t shock with that for tdyn from Figure 21 of Ivan de Ven et al 



(J2006l ) shows that tshock ~ tdyn at r ~ 27' and that tshock < tdyn at larger radii. Consequently, 



we can conclude that beyond r ~ 25-30' the energy input to the outer parts of cluster from the 
disk shocking process is significant, and it will increasingly dominate the dyna mics as the radius 



incre ases. Indeed when tshock is less than the stellar orbital timescale (~4t^y„, Ivan de Ven et al 



20061 ) the stars are unlikely to be in equilibrium with the cluster potential. It is also worth noting 
that at approximately the same cluster radius as where shock heating becomes important, the 
stellar orbital timescale exceeds the orbital period of the cluster around the center of the Galaxy. 
As a result, the outer parts of the cluster will also experience tidal heating due to the changing 
potential field of the Ga laxy as the cluster moves from its apo- to peri-Galactic distances (e.g., 
Gnedin &: Ostrikeii 119971 ) . The significance of these two effects then suggests strongly that the 
phase space structure of the outer parts of uj Cen is dominated entirely by external effects driven 
by the clu ster's location relatively close to the center of the Galaxy. The same conclusion was 
reached bv lvan de Ven et al.l (|2006l ). 



What is needed to shed further light on the situation is a full numerical simulation in which 
the dynamics of an u Cen-like system are explored as the system orbits in the potential of the 
Galaxy. Such a calculation needs to include the effects of disk-shocking and continue for a suffi- 
cient time that the quasi-e quilibrium situa,tion t hat likely applies in the outskirts of uj Cen becomes 
established. The results of iKiipper et al.l (J2010l ). for example, are suggestive in this respect. Based 
on A^-body calculat ions of model star clusters with masses of a few 10^ Mq on various orbits, 
Kiipper et al.l (J2010l ) demonstrate that tidal heating can lead to a population of "potential esca- 
pers", i.e., energetically unbound stars inside the cluster's Jacobi radius. This then results in outer 
surface density profiles that have power law slopes in the range -4 to -5. It also results in flattened 
velocity dispersion profiles that lie above the predictions of simple equilibrium models, with the 



*The value of the numerical coefficient in this equation given bv Ivan de Ven et al.l ( 20061 ) . 21 , is incorrect (G. 



van d e Venn, priv. comm. 2012). However, the error does not significantly affect the discussion in Ivan de Ven et al 
( 20061 ) . although the influence of tidal shocks at a given radius are over-estimated in that paper. 
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deviat ion commencing at about half of the Jacobi radius (JKiipper et al.ll20ld ). While iKiipper et al 
([20101) caution that their results are not readily scaleable to more massive globular clusters (recall 
the mass of w Cen is a few 10^ M©), their model calculations are at least qualitatively in agreement 
with the u Cen observations. 

More specific progr e ss in this direction is given by the results of the N-body calculation de- 
scribed in ISollima et alj (|2009l ). This is a A^ = 50,000 particle model for the cluster, i.e., central 
concentration, tidal radius and mass similar to the real cluste r, calculated for an o j Cen-like orbit in 
a three component (bulge+disk+halo) Galactic potential (see lSollima et al.l (J2009l ) for details). The 
calculations covered '^10 orbits (as against the many lO's of orbits made by the real cluster over 
a Hubble time) . The velocity dispe rsion profile of the bound remnant at the end of the simulation 
(see Fig. 11 of ISollima et al.l ()2009l )) is shown in Fig. llOi The agreement with the observations is 
excellent and it again suggests there is no need to invoke dark matter or non-Newtonian gravity to 
explain the observed velocity dispersion profile. 

In summary then, the new observations presented here confirm that the velocity dispersion 
profile of uj Centauri remains relatively fiat at ~6.5 km s~^ beyond approximately 25' from the 
cluster center. The most likely explanation of this effect is that we are seeing the consequences 
of external infiuence on the dynamics of the outer parts of the stellar system, which contain only 
a small fraction of the cluster stellar mass. Consequently, there is no requirement to invoke the 
presence of dark matter or non-standard gravitational theories to explain the observations. 
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